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Abstract 

We study the proposal that mass-varying neutrinos could provide an explanation for the 
LSND signal for — > v e oscillations. We first point out that all positive oscillation signals 
occur in matter and that three active mass-varying neutrinos are insufficient to describe all 
existing neutrino data including LSND. We then examine the possibility that a model with 
four mass-varying neutrinos (three active and one sterile) can explain the LSND effect and 
remain consistent with all other neutrino data. We find that such models with a 3 + 1 mass 
structure in the neutrino sector may explain the LSND data and a null MiniBooNE result for 
O.fO < sin 2 29 x < 0.30. Predictions of the model include a null result at Double-CHOOZ, but 
positive signals for underground reactor experiments and for — > v e oscillations in long-baseline 
experiments. 
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1 Introduction 



The LSND experiment has found evidence for — ► i> e oscillations at the 3.3a level P3E], with 
indications for v„ — * v e oscillations at lesser significance [21 El- The combination of the LSND 
data with the compelling evidence for oscillations in solar [1J H3 d, atmospheric [HI 03 El EEU 
112) . accelerator |13j . and reactor ^1] neutrino experiments cannot be adequately explained in the 
standard three- neutrino picture |15j . Extensions to models with four light neutrinos (with the extra 
neutrino being sterile) ^Hj ° r CPT violation with three neutrinos [171 118j have been proposed to 
accommodate all neutrino data. However, in both cases, recent analyses indicate that neither 
scenario provides a good description of the complete data set |191 12UI . 1 The addition of a second 
sterile neutrino improves the fit substantially 22 in some models with more than three neutrinos, 
especially if there is CP violation |2Sj. Another possible solution is to have CPT violation in a 
four-neutrino model |24j . Neutrino decay has also been proposed j25 a as a means of explaining the 
LSND data. 

Mass-varying neutrinos 26 (MaVaNs) have been discussed as a means for generating dark en- 
ergy [27j, explaining the LSND results [2S] and have been shown to improve the neutrino oscillation 
fit to solar neutrino data [22]. It is a relevant fact that for all positive oscillation signals (solar, 
atmospheric, K2K, KamLAND, and LSND) the detected neutrino travels through matter for some 
or all of its path length. In all but solar neutrino oscillations the path lengths are nearly all through 
Earth matter, while in the solar case the important matter effects occur in the sun. Therefore if 
MaVaNs exist, it is possible (although not required) that the masses and mixing angles indicated 
by the data could differ substantially from their vacuum values. 

There are three positive oscillation signals in which the neutrino path is primarily in the Earth's 
crust: KamLAND, K2K and LSND. The latest data from K2K T3] now yield an allowed region 
for the mass-squared difference which is roughly similar to that obtained for the oscillations of 
atmospheric neutrinos and is inconsistent with the mass scales indicated by KamLAND and LSND 
neutrino oscillations. The usual argument used to exclude three-neutrino models from solar, at- 
mospheric and LSND data ^5] is that there are only two independent mass-squared differences for 
three neutrinos. This argument applies in the MaVaN scenario to the combination of KamLAND, 
K2K and LSND experiments, all of which were conducted in Earth matter with similar density 
(and which therefore should be subject to similar mixings and mass scales). Thus a three-MaVaN 
model is insufficient to explain all of the data. 

In this paper we explore the possibility that oscillations of three active and one sterile MaVaN 

can explain all neutrino data including LSND and MiniBooNE. Since neutrinos in the KamLAND, 

x It has been suggested that certain approximations in some of these analyses have ignored small terms that 
may allow a better fit to the data, but a full study has not yet been made. 
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K2K and LSND experiments all pass primarily through Earth matter of approximately the same 
density (Earth crust), they should be consistent with the same set of mixing angles and mass- 
squared differences. If a set of oscillation parameters cannot be found that is consistent with these 
three experiments, then a four-MaVaN model is not possible. If such a set can be found, consistency 
with the solar, atmospheric and vacuum neutrino data must also be realized for the model to be 
viable. We study only 3 + 1 models, in which there is one mass eigenstate well-separated from the 
other three, since the constraints on 2 + 2 models are much stronger |19| I30|. 

For simplicity we examine a minimal 3 + 1 MaVaN scenario in which substantial MaVaN 
effects occur only for the v% and states and there is no vacuum mixing between active and sterile 
neutrinos, so that MaVaN effects are solely responsible for active-sterile mixing. We study the 
feasibility of such a model in describing current data. The MiniBooNE experiment |3 1 j is now 
taking data that will test the LSND oscillation parameters in the — > v e channel; we examine the 
consequences of this 3 + 1 MaVaN model for both positive and negative MiniBooNE results. We 
find that viable solutions exist if MiniBooNE sees no oscillations. Furthermore, in those solutions 
both the vacuum neutrino masses in the active sector and the MaVaN mass terms need not exhibit 
a hierarchy, and the mass scales for the oscillations of both atmospheric and LSND neutrinos 
are generated by MaVaN effects. Solar neutrino and KamLAND data are explained primarily 
by vacuum masses and mixings. Finally, we discuss the implications of this model for future 
experiments. 

2 A 3 + 1 MaVaN model 
2.1 Masses and mixings 

An element of the mass-squared matrix for MaVaNs in the vacuum eigenstate basis can be written 

as 

(M 2 )ij = {mAj - My) 2 , (1) 

where the rrii are the masses in an environment dominated by the cosmic microwave background and 
the Mij are the density-dependent mass terms generated by acceleron couplings to matter fields. 
We will assume that the heaviest neutrinos have masses of 0(0.01) eV in the present epoch, and 
that as a result of their non-negligible velocities, the neutrino overdensity in the Milky Way from 
gravitational clustering can be neglected J32] • The rrii (which we will refer to as vacuum masses) 
represent the masses of terrestrial neutrinos in laboratory experiments like those measuring tritium 
beta decay [SSI- We note that cosmological bounds on the sum of neutrino masses of 0(1) eV |34| 
are inapplicable to MaVaNs. Consequently, the usual relationship between neutrino dark matter 
and neutrinoless double beta decay [SS] is also rendered inapplicable. In Ref. [3^], it was pointed 
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out that so long as the acceleron does not couple to nonrelativistic neutrino eigenstates (which is 
the case under consideration), neutrino dark energy is stable. However, the stability of neutrino 
dark energy with the acceleron also coupled to matter has not been studied so far. 
We adopt a matter dependence of the form 

M .'iW=^(|)\ (2) 

where n e is the electron number density in units of iV^/cm 3 , are the values at some reference 
density and k parametrizes a power law dependence of the neutrino mass on density. In principle, 
Mij is expected to depend linearly on n e since the acceleron is assumed to evolve adiabatically and 
remain at the minimum of its potential. We allow k to deviate from unity to emphasize that a 
wider range of k is allowed by oscillation data. The choice of reference density is arbitrary; we will 
take it to be that of the Earth's crust, n® ~ 1.5. Implicit in the form of Eq. © is the assumption 
that the neutrino number density has a negligible effect on neutrino masses. Thus, it applies only 
in the current epoch when the cosmic neutrino background number density (0(1O~ 12 ) eV 3 ) is tiny. 
At earlier epochs, the neutrino number density is orders of magnitude larger and must be taken 
into account. For example, in the era of Big Bang Nucleosynthesis, the neutrino number density 
is C(10 30 ) eV 3 . Moreover, we have no reason to expect the M?- to be unaltered at earlier epochs 
since the acceleron-matter couplings may vary with redshift. 

For simplicity we will study a MaVaN model in which MaVaN effects occur only for and 
and there is no vacuum mixing between active and sterile neutrinos. The rationale for this choice 
is that if MaVaN effects are to be responsible for generating the LSND mass-squared difference 
and for mixing between active and sterile neutrinos in the Earth, then if they also involve the 
two lightest neutrinos it would be difficult to obtain the proper mass-squared difference for solar 
neutrinos. Then the evolution equations in the flavor basis for MaVaN oscillations in matter at the 
reference density is given by 
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where 

(m\ \ 
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M 2 = „ (4) 

(m3-M 33 ) 2 M| 4 

\ Mf 4 (m4-M 44 ) 2 / 

U is the neutrino mixing matrix that connects the flavor eigenstates with the mass eigenstates in 
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vacuum, given by 
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and Sj and Cj denote sin#j and cos #7, respectively, for j = s,a,x. The angle # s represents the 
usual mixing for solar neutrinos, 9 a the mixing for atmospheric neutrinos, and 6 X the mixing of v e 
with Vn in atmospheric and long-baseline experiments. The charged-current amplitude for v e — e 



forward scattering in matter is jSZI 

A = 2V2G F n e E u = 1.52 x l(T 7 eV 2 n e E u (MeV) 



(6) 



There is also a v — e neutral-current forward scattering contribution for all active neutrinos given 
by — Ar, where r = n n /(2n e ) and n n is the neutron number density in units of A^/cm 3 , which is 
equivalent to a +Ar term for the sterile neutrino (see Eq. [31). 

It is convenient to parametrize the mass-squared matrix in matter in terms of its eigenvalues 
M 2 and the mixing angle 9 that diagonalizes it, 
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then the mass-squared differences are SMf- = M? — Mj. If the ordinary matter term can be ignored 
(which will be approximately true for all but solar neutrinos), the neutrino mixing between the mass 
eigenstates in matter and the flavor eigenstates is 



V = U 
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where s = sin^ and c = cos 9. For solar neutrinos, where the ordinary matter terms are important, 
V cannot easily be written in closed form. Since the MaVaN terms do not affect the first and 
second generations, M\ = m\ and M2 = ui2- 

In 3 + 1 models there is one neutrino mass well-separated from the others by the LSND mass- 
squared difference (5m 2 L ), and the sterile neutrino couples strongly only to the isolated state. The 
atmospheric and solar mass-squared differences will be denoted by 5m 2 and 5m 2 , respectively. 
There are four possible mass spectra in 3 + 1 models, depending on whether the isolated state is 
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above or below the others, and whether the other three neutrino states have a normal (5m 2 > 0) or 
inverted (5m 2 < 0) mass hierarchy. We only consider the case with M4 > Mi, M%, M3 and normal 
hierarchy, which implies <5M| 3 ~ 5M\ 2 ~ 5M\ X = 5m\ > <5Mf 2 ~ SM^ = 5m 2 > 8M^ 1 = 8m 2 s . 

For now we take mi = Mi = 0, so that m\ = M 2 = 5m 2 ., M| = 5m 2 and M 2 = 5m\ (we 
discuss the case mi 7^ in Sec. 4). We will also set m^ = 0, i.e., the sterile neutrino is massless in 
vacuum. In our illustration we take sin 2 26 s = 0.80 and sin 2 29 a = 1.00, and allow 6 X to vary. For 
the remaining parameters (m^, M33, M34, M44) we have the following relations that follow from the 
diagonalization of M 2 (Eqs. ElandEJ): 

(mg - M 33 ) 2 = c 2 M 2 + s 2 M 2 , (9) 
M| 4 = s 2 Ml + c 2 M\ , (10) 
M 3 2 4 = sc(M 3 2 - M\) . (11) 

We will also assume that the vacuum masses mi and 7713 are much smaller than the MaVaN 
parameters in the Earth's crust (M33, M44 and M34), which will be justified by our numerical 
results in Sec. 3. 

In order to determine the allowed MaVaN parameters, we will use the positive oscillation results 
from KamLAND, K2K and LSND. Since these experiments were conducted in the Earth's crust, the 
mixing matrix V for all of them should be nearly the same; we will consider it as the same matrix 
for all three experiments and attempt to determine the parameters from the combined data. Since 
the baselines for these experiments are all 250 km or less, ordinary matter effects due to coherent 
forward scattering are small and the A term in Eq. © can be ignored to a good approximation. 
Therefore the mass-squared differences relevant for oscillations in these experiments are <5M 2 . and 
the mixing matrix is V in Eq. 

There are also constraints from other experiments, but only those experiments which were 
conducted primarily in Earth matter are relevant for the matrix that describes the results of Kam- 
LAND, K2K and LSND. The CHOOZ j^Sj reactor constraint on v e — ► v e oscillations at the Sm 2 
scale (L/E u ~ 250 m/MeV) does not apply since the neutrino path in the CHOOZ experiment 
was primarily in air; it is instead replaced by the weaker Palo Verde [39_ constraint at a somewhat 
smaller L/E u value. Similarly, bounds from the Bugey reactor experiment [1U] at the bm\ scale 
{L/E v ~ 25 m/MeV) must be replaced by the considerably weaker bounds from Gosgen jlj and 
Krasnoyarsk |42| . The CDHSW |43| bound on ZAj — > v„ oscillations at the 5m 2 L scale applies to the 
Earth matter the neutrino path was approximately 90% in matter in this experiment |44j . 

Once the MaVaN parameters have been determined from the Earth crust data, consistency with 
the atmospheric, solar and vacuum data can be checked. 
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2.2 Oscillation formulas 

In this section we list the oscillation probabilities in the limit that the ordinary matter effect can be 
ignored, except for solar neutrino oscillations. The relevant oscillation probabilities in the leading 
oscillation are approximately 

P{pH > ^e)LSND * 4 S 4 |[/ M3 | 2 |£/e3| 2 sin 2 A L , (12) 

i 3 (^^^)cDHSW ^ l-4s 2 |^ 3 | 2 (l-s 2 |^3| 2 )sin 2 A L , (13) 

P(^^^e)Go 8 gen - 1 " 4 S 2 |C/ e3 | 2 (l - S 2 |C/ e3 | 2 ) sin 2 A L , (14) 

where A^ is the largest of the usual oscillation arguments Aj = 5m 2 L/(4E u ) for j = L, s, a. 
The relevant oscillation probabilities at the first subleading scale (5m 2 ) are approximately 



P(^^^)K2K * l-2 S 2 |C/ M 3| 2 (l-s 2 |^3| 2 )-4c 2 |C/ M3 | 2 (l-|^3| 2 )sm 2 A a (15) 
P(v^Ve)K2K - 2 S 4 |C/ e3 | 2 |C/ At 3| 2 + 4c 2 |C/ e 3| 2 |C/ / ,3| 2 sin 2 A a , (16) 
P(^-»^)k2K * 2s 2 c 2 \U^\\ (17) 



P(v e ->• f e )paloVerde - sin 2 29 x 



-s 4 + c 2 sin 2 A a 
2 



(18) 



where we have averaged over the leading oscillation scale. We do not consider CP violation, so the 
oscillation probabilities for neutrinos and anti-neutrinos are the same. 
At the smallest scale (5m 2 ) we have 

P(Ve ^KamLAND ^ 1 " 2| U e3 \ 2 (1 - | U e3 \ 2 ) - 2 S V \U e3 \ A - 4\U el \ 2 \ U e2 1 2 sin 2 A s , (19) 

where we have averaged over the oscillations of the two higher 5m 2 scales. For solar neutrinos, for 
the large mixing angle (LMA) solution with adiabatic propagation we have 

4 

P(u e ^u e ) = J2\U ej \ 2 \V e ° j \ 2 , (20) 

3=1 

where V° is the corresponding value of V at the point in the sun where the neutrino is created 
(n e ~ 80). The fraction of solar neutrinos that oscillate to sterile neutrinos is given by 

4 

P(v e ^v s ) = J2\Uej\ 2 \V s /. (21) 
3=1 

3 Model constraints 
3.1 Fits to Earth crust data 

As discussed above, KamLAND, K2K and LSND should all be described with the 5M 2 - and mixing 
V that exist in the Earth's crust (n e ~ 1.5). KamLAND data ^1] give 5m 2 ~ 8 x 10~ 5 eV 2 with 
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oscillation amplitude 0.8 and the K2K results |13j imply 5m 2 , ~ 2.8 x 10 -3 eV 2 with maximal 
mixing (both in two-neutrino fits). In a three-neutrino model, if the mixing angle 9 X vanishes, 
then these two oscillations decouple from one another and the two-neutrino fits may be used 
directly. However, the amplitude of the LSND oscillation is (from Eq. I12j) 

sin 2 26» L = sS 2 sin 2 26^ , (22) 

so a non-zero 9 X is required to generate an LSND signal. We will examine solutions with small 
sterile mixing (sin 2 9 <C 1), so that reactor constraints at the 5m 2 scale are similar to the usual 
three-neutrino case (see Eq. 118(1 . Similarly the oscillation probabilities for KamLAND and solar 
neutrinos are approximately the same as in the three-neutrino case (see Eq. 1191 and Section 3.3). 
There are upper bounds on 9 X from three-neutrino fits |46( . sin 2 29 x ^ 0.19 (0.28) from solar and 
KamLAND data and sin 2 29 x < 0.16 (0.23) from CHOOZ, atmospheric and K2K data, where the 
bounds are at the 2a (3a) level. However, for MaVaNs traveling through the Earth's crust the 
CHOOZ bound is no longer applicable and must be replaced by the Palo Verde constraint, which 
is less restrictive by approximately a factor of two at the Sm 2 indicated by K2K and atmospheric 
neutrinos (~ 2 — 3 x 10 -3 eV 2 ). Therefore, a value of sin 2 29 x as large as ~ 0.3 would appear to be 
allowed at the 3a level with the CHOOZ constraint removed. 

Finally, it is well-known that a combination of reactor and accelerator constraints disfavor the 
standard 3 + 1 model (191 147j . In our 3 + 1 MaVaN model, the strong bounds from the Bugey 
reactor are replaced by the considerably weaker bounds (by approximately a factor of three) from 
Gosgen and Krasnoyarsk. In the region of interest {5m\ ~ 1 eV 2 ) the upper bound on the — ► 
oscillation amplitude from CDHSW is about 0.1; the corresponding bound on the v e — * v e oscillation 
amplitude from Gosgen/Krasnoyarsk is also about 0.1. From Eqs. l(13j) and ((14(1 we see that in this 
model the oscillation amplitude for Gosgen/Krasnoyarsk is smaller than that for CDHSW by a 
factor {Ues/U^l 2 = t&n 2 9 x /s 2 (for small 9), which for sin 2 29 x < 0.3 and 9 a = 7r/4 is of order 
0.15 or less. Therefore if the CDHSW bound is satisfied, then the Gosgen/Krasnoyarsk constraints 
are automatically satisified, and we need to consider only the effect of the CDHSW bound on the 
model parameters. 

In Fig. ^]we show regions in LSND amplitude sin 2 29 l and mass-squared difference 5m\ ac- 
cessible to the model and the constraints from data. The solid curves are the upper bounds on 
sin 2 29l allowed by the 90% C.L. CDHSW constraint for several values of 9 X . The intersections of 
the dashed curves with the solid curves show the values of sin 2 9 required to achieve those upper 
bounds; i.e., s = 2 sin 29 L / sin 2 29 x from Eq. (j22J). The shaded regions show the current 99% C.L. 
allowed region from a joint analysis 48] of LSND and KARMEN gHj data. For sin 2 29 x < 0.30, 
sin 2 9 > 0.23 is required to be consistent with the LSND /KARMEN 99% C.L. allowed region. How- 
ever, as we will show in the next section, such large values of 9 are not compatible with atmospheric 
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0.001 0.01 0.1 

sin 2 26 L 



Figure 1: Regions in <5m 2 -sin 2 space accessible to the 3 + 1 MaVaN model discussed in this 
paper. The solid (dashed) curves show the upper bound on the LSND amplitude sin 2 29 l from 
CDHSW for given values of sin 2 29 x (sin 2 9, the mixing between the third and fourth mass eigen- 
states due to MaVaN effects). The intersections of the dashed curves with the solid curves show 
the values of sin 2 9 and sin 2 29 x required to achieve the maximum LSND amplitudes allowed by 
CDHSW. The shaded regions show the current 99% C.L. region allowed by a joint analysis jlH] of 
LSND and KARMEN data. 

neutrino data. 

3.2 Atmospheric neutrinos 

Expressions for the oscillation probabilities for atmospheric neutrinos are similar to those for K2K 
(Eqs. H5ll7j) . except that the values for 5M 2 and V vary as the electron number density along the 
neutrino path varies (there is also an additional matter effect for v„ — * u e , v s oscillations for the 
higher-energy atmospheric neutrinos, which we ignore). We will call the varying mass and mixing 
parameters for atmospheric neutrinos M|, M 2 and 9. These quantities obey relations similar to 
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Eqs. ©-(HQ) 



(m 3 - M 33 r k ) 2 = c 2 Ml + s 2 M\ , 


(23) 


M^r 2k = s 2 Ml + c 2 M| , 


(24) 


Mf 4 r 2fc = sc(M| - Ml) . 


(25) 



with M3, M4 and # replaced by their tilde counterparts and with the MaVaN parameters M33, M34 
and M44 multiplied by the factor r k , where r is the ratio of the average matter density for a given 
path compared to the density of the Earth's crust; e.g., for a path through the center of the Earth, 
r ~ 3. For small 7713 (the solutions we are investigating) the value of the varying sterile mixing in 
the Earth, 9, is very similar to the sterile mixing in the crust, 9; this can be seen by comparing the 
expressions for 6 and 8 

?M 2 ~ ?M 2 r 2k 

tan20 = - . , tan2c? = 7 _ . (26) 

(m 3 - M 33 ) 2 - M| 4 (m 3 - M 33 r fc ) 2 - M 2 4 r 2fc V ; 

In the limit that |m 3 | <C I-M33I, it is evident that 9 ~ 9. Therefore if 9 is small then 9 will also be 
small, and sterile mixing will not upset the atmospheric neutrino fits. In this same limit the size 
of M 2 , i.e., the largest varying oscillation mass scale in the Earth, is approximately r 2k 5m 2 L . 

For 9 a = 7r/4, the amplitude for oscillation to sterile neutrinos is given approximately by s 2 c 2 
(from the tilde equivalent to Eq, I17|) : these oscillations occur at the leading mass scale (^1 eV 2 ) 
and oscillations with amplitude of order s 2 c 2 are seen for downward as well as upward neutrinos. 
Since no large suppression of downward events is observed ^0], a value of sin 2 9 (and hence also 
sin 2 9) greater than about 0.10 is disfavored. As noted in the previous section, if sin 2 29 x < 0.30, 
then the sterile mixing angle must satisfy sin 2 9 J> 0.23 to obtain a value for the LSND amplitude 
consistent with the 99% C.L. allowed region from LSND and KARMEN. Therefore, if MiniBooNE 
(in which the neutrino path is primarily in Earth crust) were to confirm the LSND/KARMEN 
99% C.L. allowed region, our 3 + 1 MaVaN model would be disfavored. 2 

The situation changes if MiniBooNE reports a null result. The lighter shaded region in Fig. [21 
shows the region that would be allowed (at 99.5% C.L.) by a combination of LSND and a null 
MiniBooNE result with 5 x 10 20 protons on target (P.O.T.) pDj. This LSND/null- MiniBooNE 
region is shifted to smaller sin 2 29 'l compared to the LSND/KARMEN region in Fig. ^ This shift 
occurs because a null MiniBooNE result would be in conflict with LSND, and a combined fit to the 
two experiments essentially results in a weighted average of the two oscillation probabilities (0.25% 
from LSND and 0% from MiniBooNE, respectively). Since LSND and MiniBooNE would be in 
conflict, no region is allowed below the 98% C.L. [HOf. As is evident from Fig. [21 this new allowed 
region could comfortably be explained by our 3 + 1 MaVaN model with both sin 2 29 x < 0.30 



2 We note that this does not necessarily rule out a more general four-neutrino MaVaN model. 
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Figure 2: Same as Fig.0 except that the lighter (darker) shaded region shows the region that would 
be allowed at 99.5% CL if MiniBooNE sees no signal after running with 5 x 10 20 (10 21 ) P.O.T. 

and sin 2 8 < 0.10, in the region near 5m? L ~ 1 eV 2 and sin 2 28l ~ 0.0003 — 0.001. An increase 
to 10 21 P.O.T. in MiniBooNE decreases the size of the combined allowed region but does not 
eliminate it (see the darker shaded region in Fig. EJ) ; there is no region below the 98.8% C.L. with 
10 21 P.O.T. [HOI. The smallest value of sin 2 28 x consistent with the LSND / null-MiniBooNE region 
is about 0.10 if sin 2 6 is not allowed to be larger than about 0.10, and sin 2 6 must be larger than 
about 0.04; see Fig. EJ 

So far we have considered only the effects of sterile mixing on atmospheric neutrino oscillations, 
which are more or less independent of the exact neutrino masses in matter, since the effects primarily 
depend on the size of 6 and not the actual mass-squared differences. However, the mass-varying 
mass-squared difference that drives the oscillations of atmospheric neutrinos (M3) will be different 
for different neutrino paths through the Earth. We need to show that the model can simultaneously 
give the correct Sm 2 for K2K in the crust and for atmospheric neutrinos. 

Although a detailed analysis of atmospheric neutrinos would be required to determine the precise 
effects of the density profiles on the allowed regions of the parameters, some semi-quantitative 
statements can be made. Since the largest oscillation signal occurs for upward events, which pass 
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through the core, as a first approximation we consider a path through the center of the Earth as 
representative of the atmospheric data; the average matter density for these events is r ~ 3 times 
the density in the Earth's crust. There are then a total of six potential observables from Eqs. Q- 
(HU) and P3>(|3K|) : Mf,Mf and 8 (from K2K and LSND), and Mf,Mf and 9 (from atmospheric 
neutrinos). 

Since the sterile mixing angle 8 for atmospheric neutrinos which pass through the core is very 
similar to 8 in the crust, its exact value is unimportant, as long as it is small. Furthermore, the 
value of the largest mass scale for atmospheric neutrinos which pass through the core, M 2 , is not 
determined by data since no such oscillations are observed. Therefore the only relevant observable 
from the upward atmospheric data is the effective mass-squared difference Mf ; we can rewrite 
Eqs. (|25 |1 -(P5 |) in terms of Mf as 

M 3 4 4 r 4fc + Mf [(m 3 - M 33 r k ) 2 + M 2 4 r 2fc ] = (m 3 - M 33 r k ) 2 M 2 4 r 2k + M 3 4 . (27) 

There are four free parameters in Eq. (j3J): m 3 , M 33 , M34 and M44. Equations (|9jl-pij) and (|2*T|) can 
be used to determine these four parameters using the oscillation data as follows: Mf = Mf = dm 2 
(assuming K2K and upward atmospheric neutrinos have the same 5m 2 ) and Mf = Sm 2 L , where 
particular values for the vacuum mixing U and sterile mixing angle in the crust 8 are also used as 
inputs. Once m 3 , M 33 , M44 and M34 are determined, the size of the LSND amplitude in Eq. (|22|) 
can be checked for consistency with the LSND result. 

Using 5m 2 = 2.0 x 10 -3 eV 2 for both K2K and atmospheric neutrinos, Table ^ shows the 
maximum value of sin 2 8 allowed by CDHSW, the corresponding LSND amplitude sin 2 2#£, and 
the MaVaN parameters m 3 , M 33 , M34 and M44 for given values of bm\ and sin 2 29 x with k = 1 
and r = 3. We see that the sin 2 28 l values are in the range allowed by the combined LSND/null- 
MiniBooNE fit (see Fig-EJ). The value of m 3 is small, about 0.006 eV, which confirms our assumption 
that \m 3 \ <C IM33I. Since the value of m<i is y 5m 2 . ~ 0.009 eV, then m 3 ~ mi, i.e., all of the 
non-zero vacuum masses have similar size. Likewise, all the MaVaN parameters M33, M34 and M44 
have similar size, of 0(1) eV. Therefore no hierarchies in either the vacuum masses of the active 
neutrinos or MaVaN parameters are required to achieve the appropriate masses and mixings. 

We note that there is some fine tuning required so that the K2K and atmospheric neutrino 
mass-squared differences are small compared to the LSND scale. This amounts to requiring \m 3 — 
M33IM44 ~ Mf 4 ; since \m 3 \ <C IM33I, this is approximately equivalent to requiring IM33M44I ~ Mf 4 . 
Since the MaVaN parameters scale similarly with density, this fine tuning is not greatly upset by 
changes in density, e.g., from the Earth's crust to its core. This is confirmed by Fig. |2J which 
shows the value of Mf for the first example in Table ^ as r is increased from 1 to 3. Since Mf 
and 9 do not vary greatly over the entire range of Earth densities, the oscillation formulas in 
Sec. 2.2 should provide a good approximation to the actual oscillation probabilities. Furthermore, 



12 



Table 1: Maximum value of sin 2 9 allowed by CDHSW, the corresponding LSND amplitude sin 2 29l, 
and the mass parameters m^, M33, M34 and M44, for several sets of values for 5m 2 L and sin 2 2^. 
The density exponent factor used in Eq. PJ) is k = 1. and the ratio of core path density to crust 
path density used in Eq. (|2?j) is r = 3. 



bm\ sin 2 29 x (sin 2 9) max 
(eV 2 ) 


sin 2 29 L 


m 3 M 33 M 34 M 44 
(eV) (eV) (eV) (eV) 


1.0 0.30 0.063 
1.0 0.20 0.061 


5.9xl0" 4 
3.7xl0~ 4 


0.0056 -0.249 0.492 0.968 
0.0057 -0.245 0.489 0.969 


0.9 0.30 0.074 
0.9 0.15 0.071 


8.2 xlO" 4 
3.8xl0~ 4 


0.0055 -0.256 0.485 0.913 
0.0056 -0.251 0.480 0.914 


0.8 0.30 0.085 
0.8 0.15 0.081 


l.lxlO" 3 
4.9xl0~ 4 


0.0056 -0.259 0.472 0.856 
0.0057 -0.252 0.467 0.858 


0.7 0.20 0.100 
0.7 0.10 0.098 


l.OxlO" 3 
4.8xl0~ 4 


0.0056 -0.262 0.458 0.794 
0.0056 -0.260 0.456 0.795 



the model should provide good fits to the atmospheric neutrino data over all zenith angles in which 
the neutrino path is mostly in the Earth. The parameters M33,M34 and M44 are also changed 
very little by changes to the input values for k (the power that determines the density dependence 
of the MaVaN parameters); r/13 increases by about 10% if k is decreased to 1/2. Therefore the 
general features of the model are insensitive to the exact density dependence of the MaVaN effects, 
to changes in the matter density as the neutrinos traverse the Earth, or to the precise value of the 
core/crust density ratio. 

The pathlength in Earth matter for atmospheric neutrinos is given by 

L m = \J B 2 cos 2 9 + sin 2 ®{2Re - e 2 ) -(R-e) cos G , (28) 

where is the zenith angle (zero for downward events), R is the Earth's radius and e is the detector 
depth (of order 1 km). For comparison, the total path length is 

L = ^R 2 cos 2 G + sin 2 Q(2Re - e 2 ) + 2R5 + 5 2 -(R-e) cos , (29) 

where 5 is the height above the Earth's surface at which the neutrino is created (of order 20 km). If 
L m is a few oscillation lengths or more, the oscillations average and the distinction between L m and 
L is inconsequential. However, when L m is significantly different from L and the oscillations do not 
average, the matter mixing angles no longer accurately describe the oscillations. For example, for 
horizontal events (G = vr/2) the total distance is less than an oscillation length for the atmospheric 
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Figure 3: Variation of M| (which determines the <5m 2 for atmospheric and long-baseline neutrinos) 
with electron number density for the first example listed in Tabled Similar behavior is exhibited 
for the other examples. 

5m? and L m /L ~ + e) — v/1/20 — 0.22. For events with near 7r/2, the difference between 

L m and L changes rapidly with 0; a detailed analysis of atmospheric data is required to determine 
the precise constraints on the MaVaN parameters. 

If 6 a vr/4, i.e., the atmospheric neutrino mixing is not maximal, then the maximum LSND 
amplitude for a given value of 9 X can be larger, as given by the formula 

sin 2 26T ax = L Wa - n/A . (30) 

L l-cos26» a v ; 

For example, for sin 2 29 a = 0.9 the maximum LSND amplitude increases by about 50%. This has 

the effect of shifting the solid curves to the right in Figs. ^ and [21 by the same factor. Current 

data still excludes the model, and the allowed parameters shift somewhat in the scenario where 

MiniBooNE sees a null result, but the general features remain the same. 

3.3 Solar neutrinos 

The matter density in the sun is much higher than in the Earth, but since in our model MaVaN 
effects only occur for 1/3 and ^4, SM^ has approximately the same value as in the standard MSW 
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scenario for all densities. The effects of sterile mixing are also small since v s couples only indirectly 
to the two lightest states. We have checked numerically that the MSW probability in Eq. (|20[) 
including all MaVaN and sterile mixing effects gives the usual result to within about a percent or 
less for the same values of 9 S , 9 a and 9 X . Also, we found that the fraction of solar neutrinos that 
oscillate to steriles (Eq. 121(1 is less than 1%, which is easily within the range allowed by current 
data 51 . Therefore the MaVaN mass terms and sterile neutrino mixing have a negligible effect on 
solar neutrino oscillations in our model. 

A non-zero 9 X improves the fit to the intermediate-energy solar neutrinos (compared to the two- 
neutrino case) , at the expense of a slightly worse fit to the low and high-energy solar neutrinos • 
As noted previously, a value for sin 2 29 x as large as ~ 0.30 is consistent with combined fits to solar 
and KamLAND data $jg\ and does not appear to be ruled out. The range 0.10 <: sin 2 29 x 0.30 
that would be consistent with LSND and a null MiniBooNE result is also consistent with the solar 
data. 

3.4 Vacuum constraints 

The constraints which apply to the experiments primarily in vacuum (Bugey and CHOOZ) must 
also be checked for consistency with the model. Since the vacuum m 2 is of order <5m 2 and the vacuum 
m\ = 0, all of the oscillations in vacuum occur with Sm 2 values of order <5m 2 , not accessible in short 
baseline experiments. Therefore both the Bugey and CHOOZ bounds are avoided in this model. 

We note that there are solutions to Eqs. (|^|)- (|11|) and (|27j) with 1714 = other than those 
listed in Tabled They have the same sign for and M33, with 7713 ~ M™, i.e., all parameters 
involving ^3 and are of 0(1) eV. In that case, there are vacuum oscillations at short baseline 
due to m 2 ~ 1 eV 2 . However, such solutions would have v e — > D e oscillations at short baseline with 
approximate amplitude sin 2 29 x which are ruled out by the Bugey reactor experiment. 

4 Discussion 

We conclude by discussing some of the main features of our MaVaN model with three active 
neutrinos and one sterile neutrino designed to explain the LSND data and a null MiniBooNE 
result: 

• The sterile mixing angle in the Earth's crust (9) is approximately the same as the sterile 
mixing angle in the Earth's core (9); these angles must obey sin 2 9 ~ sin 2 9 < 0.10 to agree 
with atmospheric neutrino oscillation data. 

• A large part of the region in Jm^-sin 2 29 l space consistent with the LSND/null- MiniBooNE 
region, can be reproduced in this MaVaN model with 0.10 < sin 2 29 x < 0.30 and 0.04 < 
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sin 2 6 < 0.10. The significant size of 6 X means that it should be detectable in proposed reactor 
experiments with expected sensitivity sin 2 29 x > 0.01 where most of the neutrino path is in 
Earth matter, such as Angra, Braidwood, Daya Bay |53| . or KASKA [54] . However, Double- 
CHOOZ [23, which should be sensitive to sin 2 20 x > 0.03, would see a null result since most 
of the neutrino path is in air (where the 5m 2 values are all of order 5m 2 ). The planned 
long-baseline experiments MINOS [SHI and ICARUS |57| . sensitive to sin 2 29 x > 0.05 at the 
90% C.L.JSS], should also see a positive signal in the —* v e appearance channel. The large 
value for 6 X might also be detectable in future measurements of solar neutrinos |52| . 

The generic features of the model are relatively insensitive to the precise density dependence 
of the MaVaN parameters. There is a certain amount of fine tuning between vacuum and 
MaVaN parameters required to achieve the relation Sm 2 <C <5m 2 in the Earth, but the mass- 
squared differences for atmospheric neutrinos are fairly stable under variations in the Earth 
density. 

The vacuum neutrino masses of the active states are given by mi = 0, m-2 = \J 5m 2 ~ 0.009 eV 
and 7?t,3 — 0.006 eV. If we allow m\ = Mi to be non-zero, it can be of order 7772 and 7773 if 
there is a small upward shift in the eigenvalues M| , Mj and M| by the non-zero value of m\. 
This shift does not make an appreciable difference to the parameters in Tabled I n this case 
there would be neither a hierarchy nor a zero value required in either the vacuum masses of 
the active neutrinos or in the MaVaN couplings. 

The values of 5m 2 measured in the K2K and atmospheric neutrino experiments do not have 
to be the same, since they are separate inputs in determining the model parameters. The 
central value in K2K is about 40% higher than for atmospheric neutrinos (with 90% C.L. 
uncertainties of order 50%). Although currently not significant, if the discrepancy between 
5m 2 values for K2K and atmospheric neutrinos persists, it could easily be accommodated in 
this model. 

It has been shown that MaVaN terms that involve only u\ and v<i can improve the fit to solar 
neutrino data 29 j. Diagonal and off-diagonal MaVaN couplings must be introduced for v\ 
and v% if there are no MaVaN terms coupling V\ or V2 with v% or 1/4, then the phenomenology 
at the <5t77 2 and <5t77 2 scales discussed in this paper is unaffected. The MaVaN couplings 
involving v\ and v<i would need to be about three orders of magnitude smaller than those for 
1/3 and 1/4. 

We have chosen to consider only MaVaN couplings that involve the third and fourth gener- 
ations, so that they do not impact the solar neutrino Sm 2 scale for Earth matter densities. 
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Introducing MaVaN terms that couple z^i or 1/2 to v% or z/4 can also induce — > v e oscillations 
at the LSND scale, but then at solar densities this mixing upsets the value of the solar 5m 2 . 
Although we have not performed an exhaustive parameter search, MaVaN terms that lead 
to the appropriate LSND amplitude and that also couple appreciably to vi and U2 at Earth 
matter densities appear to be problematic. 

• The MINOS experiment has recently found evidence for — > oscillations at a baseline 
of 735 km [S3] with a 5m 2 that is consistent with K2K and atmospheric data. Since the 
neutrino path in the MINOS experiment traverses matter of similar density as that in the 
K2K experiment, our model also explains MINOS data. 

In summary, we have presented a MaVaN model that can explain all neutrino oscillation results, 
including LSND, and a null result for — > v e oscillations in MiniBooNE. There is no hierarchy 
required in the vacuum masses of the active neutrinos, which are of 0(1CT 2 ) eV, and the density- 
dependent MaVaN parameters are all of 0(1) eV for the matter density of the Earth's crust. Active- 
sterile mixing is small and is generated solely by MaVaN effects. Due to the large value required 
for 6 X , the model predicts visible oscillation effects in underground reactor neutrino experiments 
such as Daya Bay and Braidwood, but a null result in the mostly above ground Double-CHOOZ 
experiment. Long-baseline experiments such as MINOS and ICARUS should see sizeable i/„ — * v e 
oscillations. 
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